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Search for serendipitous TNO occultation in X-rays 
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ABSTRACT 

To study the population properties of small, remote objects beyond Neptune's orbit 
in the outer solar system, of kilometer size or smaller, serendipitous occultation search 
is so far the only way. For hectometer-sized Trans- Neptunian Objects (TNOs), opti- 
cal shadows actually disappear because of diffraction. Observations at shorter wave 
lengths are needed. Here we report the effort of TNO occultation search in X-rays us- 
ing RXTE/PCA data of Sco X-l taken from June 2007 to October 2011. No definite 
TNO occultation events were found in the 334 ks data. We investigate the detection 
efficiency dependence on the TNO size to better define the sensible size range of our 
approach and suggest upper limits to the TNO size distribution in the size range from 
30 m to 300 m. A list of X-ray sources suitable for future larger facilities to observe is 
proposed. 
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1 INTRODUCTION 

Serendipitous Trans-Neptunian Object (TNO) occultation 
search is so far the only way to explore properties of 
kilometer-sized TNOs or smaller. Such sear ch has been 
conducted mainly in optical bands dBianco et al 



iBickerton. Kavelaars fc Welch] 120081; Rogues et al 



with few reported detections ( Rogues et al.l |200*6e 



ISchlichting et alj|2009l 120121) . In the X-ray band, occulta- 
tion by smaller objects may be detected because of a smaller 
Fresnel scale. Following an earlier discovery of putative oc- 
cultation events, allegedly caused by 100-meter size TNOs, 
in the 1996-2002 X-ray data of Sco X-l taken by the instru- 
ment Proportional Counter Arra y (PCA) on board Ross i 
X-ray Timing Explorer (RXTE) (jchang et all 120061 . [20071 ) . 
we conducted new RXTE/PCA observations of Sco X-l from 
June 2007 to Oct 2011 to clarify possible instrumental-effect 
contamination in that discovery. 

In this paper we report the final result of the effort in 
the search for serendipitous TNO occultation in X-rays with 
RXTE/PCA, taking into account the issue of detection effi- 
ciency. Earlier res ults based on da ta ta ken by Oct 2009 hav e 
been published in lLiu et all ([20081 ') and lChang et all l|201lf) . 
Currently such search is only feasible with RXTE/PCA ob- 
seravtion of Sco X-l, because of the large effective area of 
PCA an d because Sco X-l is the brightest X-ray source in 
the sky (|Chang et al]|200d . l2007h . Since RXTE was decom- 
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missioned in January 2012, new efforts have to await possible 
larger instruments in the future. 

TNOs larger than deca-kilometer size can be detected 
directly. Results of tho se observations can be found in 
iFuentes et all (|2010t ) and lFuentes. Trilling fc Holmanl (12011 1 ) 
for works using H ST, in Fraser. Brown fc SchwambT ( 201C ) 
for Subaru, and in lPetit et all (|201ll ) for CFHT. The TNO 
size distribution, from thousand-kilometer to decameter size 
and smaller, carries infromation of the collisional and dy- 
namical history of the early solar system. It is of essential 
importance to our understan d ing o f how the solar system 
forme d (see e. g. Chiang et al. (12007 ). Benavidez fc Bagatinl 
|2009l ). lFrase"i] |2009l ) and lKenvon fc Bromlevl (|2012| )V 



2 SUMMARY OF THE RXTE/PCA 2007-2011 
OBSERVATIONS 

RXTE had observed Sco X-l for many times, and yielded 
a large amoun t of data by 2 007 ( for details of RXTE in- 
struments, see I Jahoda et all ([20061)). Many millisecon d-dip 
events were found in those data I Chang et all [20061. but 
later it was pointed out that 'Very Large Events' (VLEs) 
are li kely the cause of those dips, rather th an TNO occulta- 
tions (jChang et al.ll2007i : [Jones et a!ll2008r i. VLEs are events 
that deposit more than 100 keV into an individual anode in 
a Proportional Counter Unit (PCU). RXTE/PCA consists 
of 5 identical PCUs. An instrument dead-time of 50 /j,s is set 
for each VLE during RXTE observations of Sco X-l. VLEs 
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Type 


number of trig] 


jered anodes 


average count rate 


A 


no 




1.73 ±0.29 


B 


all 




34.5 ± 11.0 


C 


more than one 


but not all 


44.7 ±9.97 


D 


only one 




9.37 ± 2.17 



Table 1. VLE types and their average count rates (in units of 
counts per second per PCU). The number of PCUs which are on 
during the observation varies from time to time. The total data 
employed in this analysis is 334 ksec, and is 1363 ksec-PCU when 
the number of PCU on is taken into account. 



Group associated VLE type significant less significant 







dips 


dips 


A 


Type A 


45 


180 


B 


Type B. no Type A 


2 


142 


C 


Type C, no Type A, B 


3 


35 


D 


Type D only 





6 


E 


no VLEs 


1 


10 



was discussed in IChang et alj (|201 lh based on the number 
of 'less significant' dips. In our search algorithm, we in fact 
rounded off the value of count deviations to the first deci- 
mal, so all those dips with deviation between —4.95a and 
—6.45cr were counted as 'less significant' dips. Since most of 
the dips are of 2-ms duration, we consider only the case of 
2-ms bins (more precisely speaking, each 2-ms bin in fact 
is 8/4096 s = 1.95 ms). There are 1.71 x 10 s '2-ms' bins in 
total in the 334-ks data. The expected number of 'less signif- 
icant' dips due to random fluctuation is 63.4 if a Gaussian 
distribution is assumed. The deviation distribution of the 
data is a dead-time-corrected Poisson distribution, which is 
about a factor of 2 ~ 3 smaller than a Gaussian at abou t 
— 5cr : see the discussion and Fig ure 1 in lChang et~aT1 l|201ll ). 
We therefore consider the total number of less significant 
dips to be 63.4/2.5 = 25.4. Among the 1.71 x 10 8 bins, 
9.1 x 10 7 bins are without the occurrence of any VLEs. The 
number of these 25.4 dips in Group E can be estimated to 
be 25.4 x 9.1/17.1 = 13.5. That of the other groups can 
be estimated with the corresponding VLE count rates. We 
therefore expect the numbers of 'less significant' dips to be 
0.23, 4.54, 5.88, 1.23, and 13.5 for Group A to E respec- 
tively. From Table [2] one can see that only that of Group 
E is consistent with random fluctuation. The numbers of 
'less significant' dips of other groups are clearly larger. The 
deviation is the largest for Group A, and gets smaller and 
smaller from Group A to D. The above estimate is only a 
rough one because (1) we consider only the case of 2-ms 
bins, i.e. treating all the dips as of duration 2 ms, (2) in 
our definition of dip groups, Type B VLEs can also occur 
in Group A dips, similarly for Type C and D VLEs, so the 
use of each VLE count rates for the estimation is not an 
accurate one, and (3) the factor 2.5 to account for the dif- 
ference from a Gaussian distribution is somewhat arbitrary. 
Nonetheless, we think the above estimate is good enough to 
indicate the VLE-association of those dips in Group A to D. 

Those significant dips are not due to random fluctua- 
tion. The above discussion strongly suggests that those of 
Group A to D are related to the occurence of VLEs and are 
therefore instrumental due to a process not yet fully under- 
stood. The 's ignificant' dip event in Group E, denoted as 
'Event El' in IChang et alj (|201rf ). is the only non-random 
dip event that shows no indication of any instrumental ef- 
fe ct. Its light curve w as fit with shadow diffraction patterns 
in lChang et al.l (|201 lh and the result suggests that it is prob- 
ably not due to TNO occultation either. 

In this sectio n we update the results presented in 
I Chang et al.1 (|201ll ) and have the same conclusion: (1) The 
fact that most of the 'significant' dip events are in Group 
A indicates that Group-A dips are instrumental. (2) The 
numbers of 'less significant' dip events of Group A, B, C 
and D are larger than expected from random fluctuation. It 
strongly suggests that VLEs, even for Type D, which trig- 
gers only one anode, can result in millisecond dip events. 
Only those dips without any association with VLEs can be 
considered non-instrumental. 

3 THE DETECTION EFFICIENCY 

In our earlier study to infer the TNO size distribution, a 
100% detection efficiency was assumed. To better under- 



Table 2. Number of dips in different groups. 'Significant' dips are 
those with the number of counts smaller than —6.5a below the 
average in an 8-second running window, where a is the standard 
deviation of the counts in each bin in the running window. 'Less 
significant' dips are those below —5. Oct but higher than —6.5(7. 

are produced by high-energy photons or particles. Whether 
they are really the cause of the millisecond dip events found 
earlier, however, was not conclusive because VLEs were only 
recorded in housekeeping data with a coarse time resolution 
of 125 ms in the data taken by 2007. In order to clarify the 
possible instrumental effect that may have cau s ed th e mil- 
lisecond dip events found in lChang et al.1 l|2006l . 12007). new 
observations with a newly designed data mode to record 
more detailed information of each individual VLE, such as 
the identification of the triggered PCU and the triggered an- 
ode and timing with 125 /xs accuracy, were conducted from 
June 2007 to October 2011, which produced in total 334 ks 
data good for the search for serendipitous TNO o ccultations. 

As already reported in lChang et al.l l|201ll ). if one dis- 
tinguishes VLEs into different types according to the num- 
ber of anodes that were triggered, one can find that indeed 
millisecond dip events are associated with VLEs of different 
types in a systematic manner. The average count rate for all 
VLEs altogether in the 334-ks data is 90.3 ± 18.8 counts per 
second per PCU. The rate for each type of VLEs is shown in 
Table [T] To see the association of dip events with VLEs, we 
divide dip events into groups according to their associated 
VLE types, as shown in Table [2] For more detailed descrip- 
tion of these classifications, as well as the al gorithm that 
we us ed to find and define the dip event, see IChang et al.l 
|201ll) . 

From Table [2] we can see that 45 of the 51 significant 
dips are of Group A. Type-A VLE only has a count rate of 
1.73 ± 0.29 counts per second per PCU. Apparently Type- 
A VLEs are very strong events that cause these dips, al- 
though exact details and why they are recorded with no 
anode triggered are not yet understood. The issue whether 
Type-B, -C, and -D VLEs can also produce millisecond dips 
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stand the detection efficieny of our approach, we conducted 
simulations to investigate the dependence of detection effi- 
ciency on the size and relative transverse speed of the occult- 
ing body. The detection efficiency surely also depends on the 
observed count rate, which is related to the brightness of the 
occulted background source and the instrument employed. 
In the followng we consider two different count rates. One 
is 8 x I0 4 cps, which is close to what we usually have for 
RXTE/PCA onservations of Sco X-l. The other is 2 x I0 6 
cps, which is about 20 times more than RXTE/PCA and 
could be achiev ed by future mission s such as the proposed 
LOFT project (|Mignani et al.ll201ll ). The orbital speed of 
the Earth is about 30 km/s. That at 40 AU is 5 km/s. The 
orbital speed of the RXTE spacecraft relative to the Earth is 
7.8 km/s. Sco X-f is in a direction at about 5.5° to the north 
of the ecliptic. The relative transever (to the line of sight) 
speed of an occulting body is usually largest when Sco X-l 
is in opposition (around end May) and smallest when the 
Earth moves towards or away from the direction of Sco X-l. 
The orientation of the motion of RXTE and of the occulting 
body will also affect that speed. We therefore consider three 
different speeds, i.e. 5, 15 and 25 km/s, in the simulation 
respectively. 

For a given count rate and a given transverse speed, 
we would like to know the chance of detecting an occulta- 
tion event caused by an occulting body of a certain size 
with our searching algorithm. Most of the RXTE/PCA- 
observed photons from Sco X-l are about 4 keV in en- 
ergy (0.3 nm in wavelength). The Fresnel scale, which is 
defined as yj Ad/2, where A is the wavelength and d the 
distance, is 30 m for A = 0.3 nm and d = 40 AU. One ex- 
pects to detect the shadow of an occulting body much larger 
than the Fresnel scale easily, while that of a body much 
smaller than the Fresnel scale will elude being detected due 
to diffraction. Another factor significantly affecting the ob- 
served shadow (light curve) is the impact parameter /?, usu- 
ally defined to be the shortest distance of the passing path 
to the shadow center in units of the occulting body radius. A 
fe w shadow diffraction patterns, comput ed with the recipes 
in iRoaues. Moncuauet fc Sicardvl (| 1987h and with a typical 
RXTE/PCA-observed Sco X-l photon spectrum, are shown 
in Figure [1] We note that with a point background source 
the diffraction patterns, expressed in units of the Fresnel 
scale, are all the same for a fixed ratio of the occulting body 
size to the Fresnel scale. That is, there is a size-distance de- 
generacy. Attempts to break this degeneracy are discussed 
in the next section. W e adopt the definit ion of the shadow 
boundary employed in iNihei et al.l (|2007T ). in which the re- 
lation between the shadow size f2 and the occulting body 



radius p is described as = 2(v / 3 3/ ' 2 + /9 3/;i ) 2/3 , in which 
both Q and p are in units of a Fresnel scale. Throughout 
this paper, when referring to the occulting body size we will 
use the diameter D, which is obviously equal to 2p. 

In our simulation, with a given count rate F, a given 
transverse speed v, and a given occulting body diameter D, 
we randomly pick an impact parameter j3 in the range from 
the shadow center (/? = 0) to its boundary (P — ^) with 
a uniform probability and a random epoch for each individ- 
ual implanted event in a 100-s segment. We first compute 
such a 100-s model light curve, then bin the light curve into 
0.25-ms bins, and then randomize the counts in each bin 
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distance (Fresnel scale) 



.0 = 2.0, 0=4.4 




-2 2 

distance (Fresnel scale) 



D = 20.0, 0=21.0 
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-15 -10-5 5 10 15 
distance (Fresnel scale) 



Figure 1. Computed shadow light curves for occulting bodies 
of different diameter D. All the 4 panels are for the case of a 
central-crossing event, that is, for a zero impact parameter 0. A 
typical RXTE/PCA observed Sco X-l spectrum is used for the 
di ffraction pattern c omputation, f! is the shadow width as defined 
in lNihei et al,l l|2007l ). D and Q are both in units of a Fresnel scale. 



j _ count rate 




diameter (Fresnel scale) 
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Figure 2. The detection efficiency of our algorithm in the search 
for serendipitous TNO occultations as a function of the occult- 
ing body diameter. Three different relative transverse speeds are 
considered. We assume these occulting bodies are at 40 AU away, 
that is, one Fresnel scale corresponds to 30 m. The upper panel is 
for the case of the typical RXTE/PCA count rate, and the lower 
one is for a much higher count rate, which may be achieved by 
future missions. 



to produce a simulated binned light curve. This light curve 
is then processed with exactly the same algorithm that we 
apply to real data. We repeat the procedure 10 times for 
each set of F, v and D, and record the number of times of 
detection, which we use to represent the detection efficiency. 

The result of the detection efficiency determination is 
shown in Figure [2] One can see that, with the current 
RXTE/PCA count rate, a detection efficiency larger than 
50% can be achieved for occulting bodies larger than 2 Fres- 
nel scale (diameter), considering the case of relative tran- 
severse speed being 25 km/s, while for the high count rate 
case (2 x 10 6 cps), that occulting body size can be pushed 
down to about 0.6 Fresnel scale. 



© 2012 RAS, MNRAS 000, [L]{7] 



4 Chang, Liu, and Chen 




diameter (Fresnel 



count rate — 2x10 6 cps 



speed = 5 km/s 
speed = 15 km/s E 
speed = 25 km/s | 




diameter (Fresnel scale) 



Figure 3. The difference between the detection efficiency times 
shadow width (r]Q) and the occulting body diameter (D) ex- 
pressed in the form of their ratio as a function occulting body 
diameter. For larger occulting bodies the fractional difference is 
small. For smaller ones whose detection efficiency is not yet neg- 
ligible it becomes large since the shadow width approaches a con- 
stant value when the size decreases. 



comes one order smaller when the object diameter decreases 
from 20 Fresnel scale to 2 Fresnel scale. For the current 
RXTE/PCA count rate, it actually becomes relatively neg- 
ligible for diameter smaller than 1 Fresnel scale. Nonetheless, 
-te as a function of D could be very steep, so the contribu- 
tion to the event rate from objects of diameter between 1 
and 2 Fresnel scale could be considerable. For the high count 
rate case, rj^lv becomes relatively negligible when the object 
diameter is smaller than 0.4 Fresnel scale. 

We estimate the upper limit to the TNO size distribu- 
tion at the level of setting N = 1 in Equation |T} . Although 
different models in the literature give different predictions 
for the TNO size distribution of sizes smaller than about 100 
km, a common trait among most models is a wavy shape, 
caused by the competetion of coagulation, erosion and shat- 
tering i n different size ranges, dow n to a s mall size of ab out 
0.1 km (iKenyon fe Bromlevll2012 l ) 0.5 k m |Fraserll2009l '). or 
1 km or so ( Benavidez fc Bag atin 200^) , below which the 
collisional equilibrium is achie ved. In the collisional equi - 
librium, one expects to have (|Q'Brien fc Greenberd 120031 1 
D~ q and 



dn 
dD 



7+P/3 
2+p/3 



(2) 




diameter (Fresnel scale) 



count rate = 2x10 6 cps 




speed — 5 km/s 
speed = 15 km/s 
speed = 25 km/s 



diameter (Fresnel scale) 

Figure 4. The product r/Qv, in units of the Fresnel scale times 
km/s, as a function of the occulting body diameter. This product 
is a relative contribution to the event rate estimate from TNOs 
of different size and different transverse speeds. 



For a background point sour ce, the even t rate can be 
estimated as (cf. Equation (8) in lChang et al.l (|201ll )) 

T d 2 



180 j 



(1) 



where N is the number of detected events, T the total expo- 
the differential size distribution (in units of 



sure time, 



number per unit length per square degree), r\ the detection 
efficiency, the shadow size, v the typical sky-projection 
relative speed, and d the typical distance to the TNOs. The 
product r/Q in the above equation was approximated by the 
diameter in Chang et al. (2011) The ratio of them as a func- 
tion of occulting body diameter D is shown in Figure [3] The 
product r]Qv represents a certain relative detection proba- 
bility. It is plotted in Figure [4] as a function of D for the 
three different speeds. We can see that the product r/Qv be- 



where p is the power index of the strength-size relation. 
In the so-called 'strength-scaled regime' (in contrast to 
the 'gravity-scaled regime' for larger bodies), p is negative. 
Most models give a q in the range from 3.5 to 4.0 for the 
collisional equilibrium. The distribution can be quite flat 
at the small-size end of the wavy-shape distribution. In 
iKenvon fc Bromley! l|2012f) q is about 1.0 in the range from 
0.1 km to 1 km. If we describe the TNO size distribution in 
terms of 



d;v 



d log D yd log D 



dn 



(£) 



-9 + 1 



(3) 



we can have, from Equation Q, 



dlogD 



' D 



1 NS jr_ )2 (4) 

,-9+1 rfil AD T V 180 W 



where v is assumed to be a typical speed independent of D. 

We consider two size ranges, that is, object diameter 
between 1.0 and 2.0 Fresnel scale and between 2.0 and 10.0 
Fresnel scale. Assuming a typical distance d — 40 AU and 
a typical relative sky-projection speed v — 25 km/s, with 
T = 334 ks and setting one detection as the upper limit, 
we have in the size range from 30 m to 60 m (recall that 
one Fresnel scale is 30 m in the current case and when- 
ever we use 'size' we mean the diameter) ( d ,1" n ) „ < 
4.0 x 10 11 deg~ 2 for q = 4.0 and 



;£)/ D= 3o 

'ttM < 3.3 x 

k d log D ) O =30 m 

10 11 deg~ 2 for q = 3.5. In the size range from 60 m to 
300 m we have ( 

q = 4.0 ' 

and 



dlogO/ D=60 



d log D I D = 6() i 



>' ' < 1.2 x 10 10 deg -2 



< 6.6 x 10 iu deg~ 2 for 
< 5.7 x 10 10 deg^ 2 for q = 3.5, 
for q = 1.0. Another 



v u iu B ^ / £>— go m 

factor of 1.14 should be applied to convert the estimate at 
the Sco X-l latitude to the ecliptic; see the discussion right 
after Equation (11) in Chang et al. (2011) These upper lim- 
its are plotted in Figure [5] 
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Sources 


Brightness 


(Crab) 


A 


/3 


Other names 




4U 


Al 


ASMquick 








/ITT 1 Q OK 

4U 1/oo-zo 


1.21 


0.56 


0.58 


269.56 


-1.63 


itji'zkq okh { ' V E 1 
Irll/Oo-ZOO, 


4U 1758-20 


0.63 


0.37 


0.30 


269.66 


2.91 


1H1758-205, GX9+1, Sgr X-3 


4U 1617-15 


17.95 


9.1 


11.0 


245.14 


5.73 


1H1617-155, Sco X-l 


4U 1813-14 


1.00 


0.33 


0.37 


273.24 


9.35 


1H1813-140, GX17+2, Ser X-2 


4U 1702-36 


0.79 


0.41 


0.44 


258.11 


-13.50 


1H1702-363, GX 349+2, Sco X-2 


4U 1642-45 


0.48 


0.18 


0.25 


255.31 


-23.06 


1H1642-455, GX340+0 


4U 1837+04 


0.30 


0.14 


0.12 


280.60 


28.09 


1H1837+049, Ser X-l 


4U 2142+38 


0.58 


0.47 


0.38 


345.98 


47.96 


1H2142+380, Cyg X-2 


4U 1956+35 


1.24 


0.29 


0.42 


312.86 


54.26 


1H1956+350, Cyg X-l 



Table 3. Non-extended bright X-ray sources, which are potential background targets for serendipitous TNO occultation search in 
X-rays. All the 9 listed sources are brighter than 0.1 Crab in all the 4th Uhuru (4U), HEAO 1 A-l (Al), and RXTE All Sky Monitor 
(ASMquick) catalogues. The 4U catalog was compiled with observations conducted during 1972-1973, the Al catalog was during 1977, 
and the ASMquick information is based on weekly average in August 2011. Those listed in the 2nd, 3rd and 4th columns are source 
brightness in the 4U, Al and ASMquick catalogues respectively, all in units of a Crab. Their positions are in the 5th and 6th columns, 
where A is the ecliptic longitude and the ecliptic latitude, both in degrees. This list is sorted with the absolute value of the ecliptic 
latitude 0. 



4 SUMMARY AND DISCUSSION 

We conclude the current effort of serendipitous TNO oc- 
cultation search in X-rays with upper limits to the TNO 
size distribution in the size range from 30 m to 300 m. In 
RXTE/PCA observations of Sco X-l from June 2007 to Oc- 
tober 2011, only one non-instrumental dip event was found 
in the 334-ks data. Due to the distance-size degeneracy in 
the occultation light curve, one cannot easily determine the 
occulting body size or distance. One way to go around is 
to fit the light curve with different diffraction patterns and 
find the best- fit result of the relative transverse speed, which 
can be obtained only in units of Fresnel scale per unit time. 
Then with some orbital assumptions, this speed can be com- 
pared to the real relative transvers e speed of a body at a 
certain dist ance to find the ma tch l|Chang et al.l [201ll ) . As 
discussed in IChang et all (|201lt ). based on the fitting result 
of the aforementioned approach, that non-instrumental dip 
event might be due to a TNO of 150-m size, but with a rare 
retrograde orbit, or, it might be due to an MBA of 40-m 
size, but the implied size distribution is incredibly steep, or, 
it might be due to a very nearby object of meter size moving 
at a relative speed of a few kilometers per second. The main 
reason of why there exists such a huge uncertainty in the re- 
sult is the small number of photons detected at millisecond 
time scale. 

Another possible way to break the distance-size degen- 
eracy is to relax the assumption of the background star be- 
ing a point source. The size of the X-ray emitting region 
in Sco X-l is still under debate, ranging from about 50 km 
to 50,000 km, depending on different models. Sco X-l is 
at a distance about 2.8 kpc. The correspondig angular size 
of its X-ray emitting region is 0.1 Fresnel scale at 40 AU 
considering the 50,000 km case. It is a good approximation 
to assume a point background source in such a case. On 
the other hand, if we are confident on the knowledge of the 
emitting region size and our instrument is large enough to 
provide good photon statistics, with a given angular size of 
the background source, one can detect the difference among 
diffraction shadow patterns caused by occulting bodies of 
the same size in units of Fresnel scale at different distances, 



since the projection size of the background source at differ- 
ent distances are different in units of Fresnel scale. For Sco 
X-l, if the occulting body is in the inner Oort Cloud, say, 
at 4000 AU, its X-ray emitting region will appear to be 1- 
Fresnel scale large at that distance for the 50,000-km case. 
From detailed study of the occultation light curve, one may 
be able to distinguish the projection size of the background 
source in units of Fresnel scale and therefore to determine 
the distance. Again, this requires a larger instrument to pro- 
vide good enough photon statistics. 



After 16 years of fruitful service, RXTE was decom- 
missioned in January 2012. Future projects of building 
larger instruments have been proposed, such as Athena 
20121) . AXTAR i|Rav et al.ll20lTh . and LOFT 



IIBarcons et al 



l|Mignani et al 



2011). In particular, LOFT will be 20 times 



larger than RXTE/PCA in terms of its photon collecting 
power. In Table [3] and [4] we list 9 brightest non-extended 
X-ray sources suitable for serendipitous TNO occultation 
search in the LOFT era. X-ray sources tend to be variable. 
Those 9 sources are all brighter than 0.1 Crab in all the 
4th Uhuru, HEAO 1 Al and RXTE All Sky Monitor cata- 
logues. See table captions for more details. Those 9 sources 
are located at different ecliptic latitude, providing the pos- 
siblity to study latitude dependence of TNO distribution 
in the decameter to hectometer size range. Their estimated 
LOFT count rates are also high enough to allow investiga- 
tion in two or three energy bands, which can help to identify 
diffraction features and to determine whether the discovered 
dip events are occultation in nature. One of the 9 sources 
is a High-Mass X-rar Binary (HMXB) containing a black 
hole candidate. The other 8 are all Low-Mass X-ray Bina- 
ries (LMXBs), two of which are the so-called 'Atoll sources' 
and the other 6 are 'Z sources'. They are all neutron- star 
systems showing various rich variation in their timing and 
spectral properties. Long-term monitoring of these 9 sources 
is by itself very much rewarding for the study of accretion 
physics around black holes and neutron stars, and is at the 
same time useful for pinning down the TNO size distribution 
in the decameter and hectometer size range. 
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Sources 




LOFT count rate 




Brightness Type 




band (keV) 


count rate (cps) 


source fraction 


(Crab) 


4U 1758-25 


1-2 


4.8 X 10 3 


92.0% 


1.02 Z source, Cyg X-2 like 




2-4 


1.3 x 10 5 


99.7% 






4-6 


1.2 x 10 s 


99.7% 






6-10 


8.3 x 10 4 


99.4% 






10-20 


2.1 x 10 4 


95.4% 




4U 1758-20 


1-2 


1.3 x 10 4 


97.0% 


0.48 Atoll source 




2-4 


6.6 x 10 4 


99.5% 






4-6 


5.1 x 10 4 


99.4% 






6-10 


4.1 x 10 4 


98.8% 






10-20 


1.2 x 10 4 


92% 




4U 1617-15 


1-2 


1.8 x 10 6 


100% 


12.6 Z source (Sco X-l) 




2-1 


2.5 x 10 6 


100% 






4-6 


1.1 x 10 6 


100% 






6-10 


8.0 x 10 5 


99.9% 






10-20 


2.4 x 10 5 


99.6% 




4U 1813-14 


1-2 


6.9 x 10 3 


94.4% 


0.50 Z source, Sco X-l like 




2-4 


6.8 x 10 4 


99.5% 






4-6 


5.6 x 10 4 


99.4% 






6-10 


4.1 x 10 4 


98.8% 






10-20 


1.4 x 10 4 


93.1% 




4U 1702-36 


1-2 


1.8 x 10 4 


97.8% 


0.66 Z source Sco X-l like 




2-4 


9.1 x 10 4 


99.6% 






4-6 


7.1 x 10 4 


99.6% 






6-10 


5.5 x 10 4 


99.1% 






10-20 


1.7 x 10 4 


94.3% 




4U 1642-45 


1-2 


1.3 x 10 :i 


69.1% 


0.34 Z source, Cyg X-2 like 




2-1 


3.8 x 10 4 


99.1% 






4-6 


4.4 x 10 4 


99.3% 






6-10 


3.0 x 10 4 


98.3% 






10-20 


7.3 x 10 :i 


86.8% 




4U 1837+04 


1-2 


1.8 x 10 4 


97.9% 


0.21 Atoll source 




2-1 


3.3 x 10 4 


98.9% 






4-6 


2.0 x 10 4 


98.4% 






6-10 


1.7 x 10 4 


97.1% 






10-20 


7.1 x 10 :i 


86.4% 




4U 2142+38 


1-2 


1.7 x 10 4 


97.8% 


0.44 Z source (Cyg X-2) 




2-4 


7.6 x 10 4 


99.5% 






4-6 


4.6 x 10 4 


99.3% 






6-10 


3.0 x 10 4 


98.4% 






10-20 


8.6 x 10 3 


88.8% 




4U 1956+35 


1-2 


4.6 x 10 4 


99.1% 


0.5 black hole candidate, HMXB 




2-4 


7.3 x 10 4 


99.5% 






4-6 


4.5 x 10 4 


99.3% 






6-10 


4.4 x 10 4 


98.9% 






10-20 


3.0 x 10 4 


96.8% 





Table 4. Estimated LOFT count rates for the 9 sources listed in Table 1 . These count rates are estimated with LOFT response files 
and are based on spectral mo dels reported in the literature (4U 1 758-2 5: 1 Jackson. Chur ch fc Balucinska-Churchl j2009l ); 4U 1758-20: 
Iaria et al.l ||2005|): 4U 1617-15: iBradshaw, Geldzahler fc Fomalontl 1120031 ); 4U 2142+38: iBalucinska-Church et alj [120101 ); 4U 1956+35: 
Nowak et al.l J201 ll) : the other 4 sources: ICackett et alj ||2009| ))~ Those listed count rates are the total count rate (including background). 
The brightness, in units of a Crab, in the 5th column is based on the assumed spectral model and is for the 2-10 keV band, in which a 
Crab is 2.4 X 10 _s erg cm -2 sec -1 . 
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0.01 0.10 1.00 10.00 100.00 1000.00 

diameter (km) 

Figure 5. The TNO size distribution. Plotted here is the differ- 
ential sky surface density at the ecliptic per decade of size in units 
of number per square degree. The upper limits in the size range 
of 30-60 m and of 60-300 m, denoted with downward arrows, are 
derived from our non-detection in the 334-ks RXTE/PCA data 
of Sco X-l, and is set at the level of one detection in 334 ks. The 
cases of g equal to 3.5 and 4.0 ( dl „" D ~ D-1+ 1 ) are considered 
in the two size ranges, of which the derived upper limits are very 
close to each other as plotted in this figure. The case of q equal to 
1.0 is also considered in the size range of 60-300 m. The downward 
arrows are plotted at the size-range boundary for q equal to 3.5 
and 4.0, and at the middle of the size range for q equal to 1.0. The 
asterisk symbol at 0.5 km is bas ed on the reported HST /FGS de- 
tection of an occultatio n event {Schlichtine ct al. (20o3); se e also 
ISchlichting et al.l <2012ln . The solid curve is the double-power-law 
distrib ution of large TNOs l|Fuentes et al 1l2010h : see lChang et al] 
(I201lh for the assumptions made to convert the luminosity func- 
tion to the size distribution. The dashed line is a direct extrap- 
olation from the large size end of the double-power-law towards 
smaller size and the dotted line is a power law anchoring on the 
double-power-law at 90 km with a power index of -3.0. We note 
that the power indices in this figure, explicitly printed close to the 
corresponding power-law lines, are for the differential size distri- 
bution per decade of size, and are therefore the same as that of 
the cumulative size distribution. The double-power-law and the 
HST/FGS result in fact indicate a wavy shape for the TNO size 
distribution, similar to that of the main-belt asteroids. 
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